The many-electron ground-state determines uniquely the potential in 
Spin-Density-Functional Theory for non-collinear magnetism 
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Since Spin Density Functional Theory was first proposed, but also recently, examples were con- 
structed to show that a spin-potential may share its ground state with other spin-potentials. In fact, 
for collinear magnetic fields and systems with fixed magnetization, the mapping between potentials 
and ground states is invertible, provided the magnetization is not saturated and that spin-potentials 
are determined within a spin-constant. We complete the proof that the mapping is invertible also 
for non-collincar magnetic fields and systems with more than one electron. We then discuss the 
non-collinear exchange and correlation energy functional and suggest improvements. 
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5h 



c 

o 

> 

ON 

o> 



in 
o 
-i— > 



T3 

c 

o 
o 



X 



The remarkable success of Density Functional Theory 
(DFT) [HUSH, as a tool for electronic structure cal- 
culations in solid state physics and quantum chemistry, 
is owed to its efficiency but also to the soundness of the 
underlying principles. In fact the quest to refine the foun- 
dations of the theory often results naturally in better un- 
derstanding and improving the accuracy of the approx- 
imation. The aim of this work is to present such an 
example. 

So far, it appeared the sound basis of ground state (gs) 
DFT is not shared by gs Spin Density Functional Theory 
(SDFT) [2, Oj LJ LzI- M, III , commonly used when the elec- 
tronic system of interest lies in a weak external magnetic 
field, or even when there is no external magnetic field 
but the electrons in the system do not form closed shells. 
Although formally SDFT is a straightforward extension 
of the original theory, the analogue of DFT's Hohenberg 
Kohn (HK) theorem 0, establishing the 1-1 correspon- 
dence between the set of spin-potentials, {V}, and the set 
of gs spin-densities, {g}, has not been established yet for 
non-collinear magnetic fields. By the spin-potential we 
mean the potential and magnetic field, V = (V; B), while 
the spin-density represents the charge and magnetization 
densities g — (p; m). 

Traditionally, invertibility of the mapping between {V} 
and {g} is proven in two steps. Firstly showing that dif- 
ferent spin-potentials have different ground states, i.e., 
the mapping between the set of spin-potentials {V} and 
the set of ground states {^} is invertible. Secondly, 
showing that different ground states, arising from dif- 
ferent spin-potentials, have different densities, i.e., that 
the mapping between {\f r } and {g} is also 1-1. 

It is straightforward to show that the second map- 
ping is invertible. In the publication where von Barth 
and Hedin 5] proposed SDFT, they also constructed, for 
the single-electron gs ^ of any spin-potential, a whole 
class of different (non-collinear) spin-potentials that ad- 
mitted \& as an eigenstate, and under rather mild con- 
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ditions as their gs (see however |io|). concluding that 
the first step in the proof of the HK theorem could not 
hold in general. Recently, interest revived on the same 
issue El El El El El E3 In the case of collinear 
magnetic fields, Capelle and Vignale E3 g ave several 
examples where the spin-potential was not determined 
uniquely by the spin-density. Eschrig and Pickett E3 
studied collinear and non-collinear magnetic fields. They 
distinguished two cases, when the gs is a pure-spin (p-s) 
state or an impure-spin (imp-s) state. Pure-spin states 
are states which through a local rotation of the spin- 
coordinates may be transformed to have a definite num- 
ber of spin-up and down electrons. Impure-spin states 
are states which cannot be transformed to have a def- 
inite number of spin-up and down electrons under any 
local spin-rotation. It was found that when the gs of the 
Hamiltonian TL is an imp-s state, there is no ambiguity 
(apart from an overall constant) in determining the spin- 
potential in TL by the ground state. However they also 
discovered p-s states are eigenstates of an operator which, 
in the laboratory unrotated spin-space, corresponds to a 
non-collinear magnetic field with constant magnitude and 
concluded the spin-potential in TL cannot be determined 
uniquely by the ground state. 

If these findings did hold in general, the mcaningfulncss 
of the Kohn Sham (KS) scheme would be questioned, 
since in the latter, the non-interacting KS spin-potential 
is not known and must be determined together with the 
KS state that yields the gs spin-density. 

Recent work gives hope that the mapping is invert- 
ible. Numerical investigations by Ullrich [l|j on lattices 
give evidence that in the continuum limit and for non- 
collinear magnetic fields, non-uniqueness is very rare if 
not completely absent. Further, the study of finite tem- 
perature SDFT by Argaman and Makov [l^] , in the limit 
of vanishing temperature (provided the limit is not sin- 
gular) shows that the mapping between V and g is in- 
vertible. 

In the first part of the paper we resolve the issue prov- 
ing conclusively that the mapping is indeed invertible. 

There is no ambiguity regarding the invertibility of 
the mapping between potentials and ground states in 
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SDFT for collinear fields. It should be noted beforehand 
that for collinear magnetic fields the eigenstates of TL 
are p-s states, since the spin-up and spin-down particle- 
number operators, Nf, iVj,, independently, commute with 
TL. Capelle and Vignale 0] discovered two cases of free- 
dom in determining the spin-potential: 

First, as in DFT where the trivial density, p(r) = 0, 
cannot determine the potential, V(r), in the same way 
here a trivial component of the spin-density cannot deter- 
mine the corresponding component of the spin-potential, 
i.e., for a system of electrons that are all spin-polarized 
in one direction, say, P|(r) > 0, Pi(r) = 0, the spin- 
component of the potential in the other direction, Vj,(r), 
is completely arbitrary (as long as the system remains 
fully polarized) . This kind of indeterminacy cannot have 
any practical consequence as there are no electrons to 
experience the arbitrariness. 

Second, together with the freedom of a constant shift 
in the potential (because iVj + iVj, commutes with TL), 
there is an additional freedom of a constant shift in the 
magnetic field, B$ appearing because the magnetization 
operator N-\ — N±, commutes with TL. According to 
Capelle and Vignale, this is an example of a systematic 
non-uniqueness arising from a conserved quantity (as op- 
posed to an accidental non-uniqueness like the example 
of the perfectly spin-polarized system.) The formulation 
of SDFT can be corrected to respect this freedom by 
constraining, at the functional differentiations, both the 
charge and the magnetization densities to integrate sepa- 
rately to fixed number of particles 01 ■ Different KS cal- 
culations must be performed for systems with fixed mag- 
netizations — (fixed spin moment calculations), 
and the gs i s gi ven by the system with the lowest KS 
total energy |14j. 

Still, the challenging question in collinear SDFT, in 
order to establish the 1 — 1 mapping is whether this is all 
the freedom available, or whether there could be other 
kinds of non- uniqueness present. It was shown that for 
continuous spin-potentials this does not happen 0, 0] . 

The general case of non-collinear SDFT is not very dif- 
ferent. We focus first on von Barth and Hedin's example. 
For any single-electron state VE' they constructed a single- 
particle spin-potential V" which acting on W gives zero, 



V' * = 0. 



(1) 



Hence, if is the ground state of a spin-Hamiltonian TL, 
then \1/ will be an eigenstate of TL + A V". A whole class 
of V" was found and it was argued that small enough 
A must exist for which \P is the gs of TL + XV'. Von 
Barth and Hedin did not explain how they constructed 
V', but it is easy to do so: any single-particle state 'J 
with spinor wave- function (^j^j) defines its own (local) 

orientation in spin-space, call it spin-up, 1 1) = (^(^j) 

/x/IV'tC 1 ')! 2 + Wi ( r )l 2 - The spin-down orientation is also 

defined: ||) = ( « ) / VW^¥+¥l^W ■ Then, * 



is an eigenstate of the single-particle spin-down number 
operator \\) (J, | with eigenvalue zero and for any V'(r) the 
spin-potential V = V(r) |J.)(J. | maps "J/ to zero. It is evi- 
dent this example is a generalisation to the non-collinear 
case of the accidental non-uniqueness for systems fully 
saturated in one spin-direction, leaving the potential in 
the other spin-direction undetermined |l8j |. To extend 
von Barth and Hedin's example to many-electron sys- 
tems, one would need to construct fully spin-polarized 
p-s states and therefore to obtain a local rotation in spin- 
space that reduces TL to spin-diagonal form. 

We point out that the discovery by Eschrig and Pick- 
ett relies exactly on the assumption that this is possi- 
ble. In fact criticism on the reasoning by Eschrig and 
Pickett (Footnote 18 in 0|) misses the point, as it is re- 
stricted to Hamiltonians that cannot be transformed to 
spin-diagonal form. 

So, given the proof ^3 that imp-s states do not lead 
to any kind of non-uniqueness, the question about the 
invertibility of the mapping between spin-potentials and 
ground-states for non-collinear spin Hamiltonians boils 
down to the existence or not of a local transformation 
that spin-diagonalizes TL. We consequently prove that a 
Hamiltonian, 



H = T + XV CC + V + B . 



(2) 



cannot be transformed to spin-diagonal form, unless the 
magnetic field is collinear. 

T, V cc in are the kinetic energy and electron re- 
pulsion operators, A = 1 ( ) for an interacting (non- 
interacting) system and V, B are the operators for the 
potential and the magnetic field. We shall focus on 



T=~ E /^t (r)V 2 Vv(r) 
2m ' J 



and 



B 



dr rh(r) • B(r) 



(3) 



(4) 



where, V>|(r), ^v(r) are creation and annihilation fermion 
field operators, m(r) is the magnetization density oper- 
ator 

m(r)=- M0 (^(r)^(r)) CT (|jW) (5) 

and er is the vector of 2 x 2 Pauli spin-matrices, er = 
(a Xl (7y,a z ). 

We introduce the general unitary transformation that 
rotates locally at every point in real space the spin de- 
grees of freedom 

/ e «M cog £,,(,.) -e^Wsincj(r) \ , . 
{r) ~ \ e-'^Wsin^r) e-' w ^ cosLu{r) J { ) 

6>(r), </>(r), uj(r) are real functions. We have U(r) l)T(r) = 
Ut(r) U(r) = I. Rotated in spin-space second quantized 
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fermion field operators can be denned: 

(*}(r)*|(r)) = (^(r)^J(r))ut(r), 

Hr)) = u(r) (£W 



(7) 
(8) 



Is there a U that diagonalizes Tt in spin-space? 

A rotation in spin-space may transform locally the 
magnetic field operator B to diagonal form. V and V ee 
are independent of spin and do not change form. How- 
ever, the kinetic energy operator T in is not invariant 
under the gauge-transformation JBJ and the transformed 
operator represents the motion of a particle in a gauge 
vector potential A(r) in spin-space (omit obvious r de- 
pendence): 

f = It I dr (*} fct) HV + A] 2 , (9) 



2m 



where A(r) = -i U(r)VU t (r). 

We further consider the rotated spin-up, down num- 
ber operators Nyj T = J dr$J(r)w T (r), r =t, |, and the 
magnetization operator in rotated spin-space: 



ut 



uj 



(10) 



Rotating back to the original spin space we see that b is 
the operator of a magnetic field b(r)//iQ ®: 



<ir rii(r) 



Mr) 

Mo 



(11) 



b(r) is in general non-collinear, b x (r) = — sin2w(r) 
cos (6(t) - 4>{r)) , b y (r) = -sin2w(r) sin (0(r) - 0(r)), 
& z (r) = cos2w(r), and has unit magnitude b 2 (r) = 1. 

We obtain tautologically the argument of Eschrig and 
Pickett. A p-s state (in rotated spin-space) is by defini- 
tion an eigenstate of the rotated magnetization b, which 
in the unrotated spin-space represents a non-collinear 
magnetic field with constant magnitude. Of course, 
the question is whether p-s states exist, or equivalently 
whether an operator b exists that commutes with Ti. |17| . 
In that case, TL could be brought to spin-diagonal form. 
Indeed in Ref. [l2T |. it is explicitly assumed that such b 
exists (b is denoted U , defined by Eq. 15 in P3)- 

To answer the question we need the commutator [H, b\. 
We have [f co , b] = [V, b] = and 



\TM 



-ih 



a—x,y,z 



dr j a (r) • Vb a (r) (12) 



[B, b] = -Ai / dr m(r) 



B(r) x b(r) 



Mo 



(13) 



where, 
ja(r) 



ih 

2m 



V ($(r)$(r) 



V> T (r) N 



($(r)#(r)) o a 



, a = x, y, z 



(14) 



The two commutators correspond to operators for dif- 
ferent physical quantities, i [T, b] describes the interac- 
tion of currents with vector potentials and i [B, b] the en- 
ergy of a magnetic moment in a magnetic field and hence 
they cannot cancel each other. In order that [H,b] = 0, 
each commutator must vanish separately. From (|12|) we 
have V& Xi!/)Z (r) = 0, i.e., b is independent of r. From 
(|13|l we have B(r) x b = and B(r) must be collinear, 
as it has to be parallel to b which does not depend on r. 

Hence, in the non-collinear case, the Hamiltonian can- 
not be spin-diagonalised by a rotation in spin-space and 
for many-electron systems the mapping between spin- 
potentials and ground states is invertible 01 ■ 

Approximate collinear functionals (in the Local- 
Spin-Density Approximation, LSDA and Generalized- 
Gradient Approximation, GGA) have been employed to 
simulate the non-collinear exchange and correlation func- 
tional E xc , Hid: The magnetization density is ro- 
tated to locally spin-diagonal form: m(r) ■ cr — ► M z (r) a z 
= U(r)m(r) • erU^(r), then, the locally diagonal spin- 
magnetization M z (r), is viewed momentarily as if it were 
the magnetization density of a collinear system. The 
corresponding (exchange and correlation) collinear mag- 
netic field B z (r) is found by taking the functional deriva- 
tive of the collinear exchange and correlation functional 
E° xc : B z (r) = 5E° c [p,M z ]/6M z (r). Then, B z (r) (seen as 
the magnetic field that would have resulted by properly 
taking the functional derivative SE xc [p, M z ]/SM z (r) of 
the non-collinear functional) is rotated back, B z (r) a z — » 
B(r)-er = U^(r) B z (r) a z U(r), to obtain the non-collinear 
magnetic field B(r). Finally Pauli-like KS equations 
are solved using B(r). Perturbative corrections (spin- 
stiffness) for the non-collinear exchange energy have been 
proposed 0. The method has given access, successfully, 
to a wide range of systems exhibiting non-collinear mag- 
netis m j^j, E| , notably the spiral spin-density- wave gs of 

Next, we analyse the consequences of using the 
collinear rather than the non-collinear exchange and cor- 
relation energy for the functional derivative. The differ- 
ence is the absence of the vector potential in E® c [p, M z \. 
We shall investigate the effect this omission for weak A. 

Following Korenman, Murray and Prange 18], the 
transformed kinetic energy operator can be expressed 
as a sum of three terms, T = Tq + Ta + Va, where Tq has 
the familiar form in terms of the rotated field opera- 
tors (|7I8|) . Ta is linear in A and describes the interaction 
of A with spin-currents, Ta = T/\ x + T& y + T& z , with 



Ta 



dr J„(r) ■ A Q (r) , a = x,y,z (15) 



The spin-current operators J Q (r) are given by 
(|15fl after replacing the unrotated field oper- 
ators by the rotated ones (|7I8|I and A- X (y) — 
sin(0 + 4>)Vw - (+) cos (6» + (f>) sinw coscj V(9-<j>), 
A z = - cos 2 lu V6> - sin 2 uo V(j). 

Va contains the square of A and has the form of a 
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potential energy operator: 

Va = J drp(r)V A (r) (16) 

where, K A (r) = {h 2 /2m) (A 2 x {r) + A 2 y (r) + A 2 z (r)) and 

p(r) is the density operator J2 T ^>J(r)Vv(r). 

The non-interacting kinetic energy functional 
T s [p,m] = min$^p im ($|T|<i>) is written in rotated 
spin-space, T s [p, m] = ($ p , m |f + f A + K A |$ P ,m} , where 
®p,m is the KS state (imp-s). 

We note M z (r) represents (in rotated spin-space) the 
same magnetization density as m(r). We further define: 

T°[p,M z ] = mm (*|T |4) - {* Pl M.\f \$ PlM ,) (17) 

<P — >p,M z 

T®{p,M z ] would be the collinear non- interacting kinetic 
energy functional if M z were a collinear magnetization 
density in the laboratory spin-space. The minimizing 
(KS) state $ p m z is p-s. For weak A we can approximate 
T s [p, m] by T s ° [p, M z ] plus a correction: 

T s [p,m]=T s VA^] (18) 

+ ($p,m, I^A. |* P ,M e ) + /" dr p(r) K A (r) + AT, [p, M z ] 

A similar analysis can be carried out for the internal 
energy functional, F\p,m\ = min*- >p , m (^ r |T + K 00 ^}. 
Replacing f by T + V cc above, we obtain F°[p,M z ], 
AF[p, M z ] in place of T°[p, M z ], AT s [p, M z ] and the min- 
imizing states * p , m , * P ,Af, in place of $ p , m , $ p ,m*- 

The non-collinear exchange and correlation functional, 
E xc [p, m] = F[p, m] - T s [p, m] - (1/2) fp(r)p(r')/\r - r'| , 
becomes 

E xc [p,m] = E° c [p,M z ] + (y p , M jT A jy p> M x ) (19) 
-($p,m« |T A , |$ P>M J + A^ c [p, M,] 



where, A£ ac |>,M,] = AF[p,M z ] - AT s [p,M z ] and 
^ c [ A M z ]=^b,M z ]-T s [ P) Mj-(l/2)/p(r)p(r')/|r- 
r'|. When M z represents a collinear magnetization den- 
sity in the laboratory spin-space, then E xc [p, M z ] is the 
collinear exchange and correlation energy functional. 

In the non-collinear LSDA/GGA method 0,0], instead 
of the non-collinear E xc [p, m], one deals with E xc [p, M z ]. 
From Eq. ^1 there are corrections, (^ p ,m^ \T a ^ \ 1 i'p,M z ) — 
($p,mJ t k\®p,m*) (first-oder) and AE xc [p,M z ] (second 
and higher). These corrections are given in terms of 
p-s states <I>p : M z , ^p,m z and it is possible to approxi- 
mate them perturbatively, using the spin-polarized uni- 
form electron gas model, in order to obtain a fully non- 
collinear LSDA exchange and correlation functional. 

Recently, the formalism was developed to treat ex- 
change exactly for non-collinear magnets [2^]. Both 
approaches have advantages (simplicity of calculation, 
access to larger systems, some account of correlation 
for non-collinear LSDA/GGA, compared with accurate 
treatment of non-collinearity and elimination of self- 
interaction for exact-exchange). 

In conclusion, we have analysed the non-collincar ex- 
change and correlation energy functional. It is intrigu- 
ing to derive and observe the effect of corrections we 
have suggested and to compare results with the exact- 
exchange scheme. We remark that any KS approach, 
crucially relies on the unique determination of the KS 
potential by the KS state that this study ensures. 
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Levy for discussions and encouragement at an early stage. 
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